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Abstract. Using McMullen's Hausdorff dimension algorithm, we 
study numerically the dimension of the limit set of groups gener- 
ated by reflections along three geodesies on the hyperbolic plane. 
Varying these geodesies, we found four minima in the two-dimen- 
sional parameter space, leading to a rigorous result why this must 
be so. Extending the algorithm to compute the limit measure and 
its moments, we study orthogonal polynomials on the unit circle 
associated with this measure. Several numerical observations on 
certain coefficients related to these moments and on the zeros of 
the polynomials are discussed. 



1. Introduction 

Curtis McMullen introduced in [12] a very efficient algorithm for 
the computation of the Hausdorff dimension of the hmit set Aoo of 
general conformal dynamical systems. Taking as our dynamical system 
a group r = ... generated by reflections along k + 1 geodesies in 
the Poincare unit disk D, we study the limit set Aoo C = dV) and 
the (unique) limit measure /x supported on it. In this article, we study 
two spectral aspects of the group, its limit set and its limit measure. 

The ffist spectral aspect arises from the intimate connection between 
dim// (Aoo) on the one hand, and, on the other, the bottom of the spec- 
trum of the Laplacian in the infinite-area hyperbolic surface S, defined 
as the double cover of the quotient D/F |3l Ch. 14;|25]. In particular, 
we study the Hausdorff dimension D{do,di,d2) of a group Tg^^^g^^^g^ pa- 
rameterised by three distances do, di and ^2 with do + di + d2 = d. 
Despite our initial expectation that the global minimum of D is at- 
tained only in the symmetric configuration do = di = d2 (centre of 
triangle in Fig. [1]), numerical computations gave us three other points 
(open circles), where this global minimum is attained. This can be 
explained via a particular symmetry in the distance parameters, which 
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we subsequently proved rigorously in Prop. 14.21 Moreover, our numer- 
ical results suggest that D >l/2 whenever (ij = for some i and that 
the global maxima of D is attained when either di = d or di = dj = d/2 
(filled circles in Fig. [1]). The interplay between the group Tg^ g-^^ g^ and 
the associated hyperbolic surface S is crucial to arrive at the spectral 
results. (See also [UlII] for further studies of Laplace eigenvalues on 
families of hyperbolic surfaces.) 



B 




A 



Figure 1. Four minima (open circles) and 6 maxima 
(filled circles) of D{do,di,d2). Barycentric coordinates 
subject to do + di + d2 = d; see §4.21 for details. 

The second spectral aspect arises in connection with a family of uni- 
tary matrices called CMV matrices and orthogonal polynomials on the 
unit circle (OPUCs); see [T8l[T9] for an encyclopaedic reference. As 
discovered by Cantero, Moral and Velazquez [5] , the multiplication op- 
erator f{z) I— >■ zf{z) is represented in a suitable basis for L'^{S^]fi) 
by a semi-infinite pentadiagonal matrix bearing their names which can 
be regarded as the unitary analogue of a Jacobi matrix or a discrete 
Schrodinger operator. The characteristic polynomials of the k x k sim- 
ple truncation of this CMV matrix is a polynomial ^k{z) that forms an 
orthogonal set in L^(5'^; /i). As do their real counterparts, ^k{z) satisfy 
a recurrence relation fl3T|) . whose complex coefficients ak are known as 
Verblunsky 's coefficients. 

One example considered by McMullen in [12] is the case of three 
symmetric geodesies go, 91,92 in D, dubbed "symmetric pair of pants". 
As the opening angle 6 of each geodesic varies from to 27r/3, the 
limit measure fig describes a continuous transition from an atomic mea- 
sure supported at three points to the absolutely continuous Lebesgue 
measure on S^. Extending McMuUen's algorithm, we computed (ap- 
proximations to) the measure /le and its associated moments. These 
are then used to study the zeros of certain "paraorthogonal" polyno- 
mials $fc(z;/3), which are the k x k unitary truncations of the CMV 
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matrix depending on a parameter /3. Among our numerical observa- 
tions, we found that the Verblunsky coefficients are all negative and 
are monotonic in 9 [cf. Fig. [Sj. We also observed that the zeros of 
the paraorthogonal tend to cluster together near the gaps of 

supp/ig and they are in some sense monotone in (3 [cf. Fig. [TO] . 

Finally, we note that a different algorithm for the calculation of the 
Hausdorff dimension was introduced in f^. Baragar [1] also described 
an algorithm for the calculation of the Hausdorff dimension of three 
geodesies go, 91,92 with two distances equal to zero, for which Mc- 
MuUen's algorithm is not applicable. 

The rest of this paper is structured as follows. For the reader's con- 
venience, we first recall basic facts about hyperbolic reflection groups 
in §2.11 and describe McMuUen's symmetric example in §2.21 Descrip- 
tions of McMullen's Hausdorff dimension algorithm and its extension 
follow in Section |3l In Section HJ we prove a symmetry property of the 
Hausdorff dimension D{dQ, ^1,(^2) and present our numerical computa- 
tion of D. Section [5] describes our numerical observations related to 
OPUCs and CMV matrices, after a brief review of the known theory. 

Acknowledgements. We would like to thank Patrick Dorey, John 
Parker, Mark Pollicott and Scott Thompson for useful comments and 
suggestions. NP is grateful for the hospitality of Williams College. KG 
was supported by a Nuffield Undergraduate Research Bursary. 

2. Hyperbolic reflection groups 

2.1. Limit sets and limit measures. Let B) = {z E C : \z\ < 1} 

denote the Poincare unit disk with its hyperbolic distance function 
d{-,-) [2^, §7.2]. Geodesies are circular Euclidean arcs, meeting the 
boundary = perpendicularly. Let Cg denote the Euclidean cir- 
cle representing the geodesic 9, i.e. 9 = 3 n Cg. Then the hyperbohc 
reflection in g agrees with the restriction (to the unit disk) of the Eu- 
clidean reflection in the circle Cg. 

Let 9o, . . . , gkhe k + 1 geodesies with corresponding Euclidean circles 
Co, . . . ,Ck and closed disks Dq, . . . ,0^ C C such that Ci = dDi. We 
assume that the geodesies are not nested, i.e. the disks Dq, . . . ,Dk 
are pairwise disjoint. To this setting, we associate the discrete group 
r = Tg^^^,,,^g^ generated by the reflections pi in the geodesies gi. V acts 
on D by hyperbolic isometries. This action extends to the boundary 
5*^, and the limit set of P is given by Aoo(r) = P ■ p fl 5*^ (which is 
independent of p G D). If we have at least three geodesies, the limit 
set Aoo(P) C has a Gantor-like structure. This can be seen as 
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follows: We refer to the disks Df^ = Di as the primitive cells or cells 
of generation zero. Let pi be the Euclidean reflection in Ci. The cells 
of generation one are given by Ps{Df'^) with s ^ I. Every primitive cell 

nf^ contains k cells of generation one. The cells of generation two are 
obtained by reflecting all cells of generation one in all primitive circles 
in which they are not contained. Repeating this operation, we obtain 
successive generations of cells and a nested structure, where every cell 
of generation n > 1 has a unique parent (of generation n — 1) and 
precisely k children (of generation n + 1). If v4„ denotes the union of 
all cells of generation n, we recover Aoo(r) C as the intersection 

A F-invariant conformal density of dimension 5 > is a measure 
/i supported on Aoo(r) satisfying, for all 7 G F and all continuous 
functions / : 5*^ — t- C, 

f{z)dp= [ fo^{z)\^'{z)\'dp, (1) 

where 7' denotes the ordinary derivative of 7 : S*^ — )■ S""^ with respect 
to the angle metric in radians. This means that p behaves like a 6- 
dimensional Hausdorff measure. 

Many properties of conformal densities can be derived from the fol- 
lowing explicit construction, due to Patterson [15] and Sulhvan [2^12^ . 
Let Fq C F be the index-2 subgroup of orientation-preserving isome- 
tries (compositions of an even number of reflections) and let gs{x,y) = 
'^'y(zro^~^'^^^'"'^^ ■ Then there exists a critical exponent 5{T), indepen- 
dent of x,y G D, such that gs{x,y) is convergent for s > 6(T) and 
divergent for s < S{T). Since Fq is geometrically finite, gs{x,y) di- 
verges also at s = 5(F) [HI Thm. 9.31]. Let 

where 6z denotes the Dirac measure at z. By Helly's Theorem, there 
are weak limits of for certain sequences Sj \ 6(T) and, due to the 
divergence at the critical exponent, each such weak limit is supported 
on the boundary S^. For the following results we refer to fl^, Ch. 9] or 
to the original articles by Patterson and Sullivan: 

Theorem 2.1. Let Fq be the index-two subgroup o/F = Fgi,...,^^, Ps cls 
defined in ([2]), and (5(F) be the critical exponent. Then 

(i) every weak limit of ps is supported on Aoo(F); 

(ii) every weak limit at x = satisfies the transformation rule ([1]) 
and is therefore a T-invariant conformal density; 
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(iii) for X = y = 0, all weak limits are probability measures; 

(iv) there is only one conformal density fi of dimension S{r), up to 
scaling; 

(v) the conformal density fi of dimension 6{r) has no point masses; 

(vi) the critical exponent coincides with the HausdorfJ dimension of 

Aoo(r). 

Henceforth, we refer to the unique conformal probabihty density /xr 
of dimension 6{T) as the limit measure of the hyperbohc reflection group 

r. 




Figure 2. The geodesies for SPP with ^ = 110°. The 
primitive [red in online version] geodesies are of the same 
size, as are generation-1 [blue] ones, but not higher-gen- 
eration [green] ones. Note the 6-fold symmetry. 

2.2. McMullen's "Symmetric Pair of Pants". We now introduce 
the "Symmetric Pairs of Pants" from [12]. This example, henceforth 
SPP, will be important in Section [5] below. For 9 G (0,27r/3), let 
fi'O) fi'i) 92 be three symmetrically placed geodesies with end points Zj e^^^/^ 
for j = 0,1,2, where zj = Ce''^^ ^^"^ Ce = e^'^'/^ [the largest /red 
geodesies in Fig. |2]. For brevity, we denote the corresponding re- 
flection group by F^, and the corresponding limit set and limit mea- 
sure by Aoo{0) and fig, respectively. It is clear from the construc- 
tion that Tg respects the 6-fold symmetry generated hj z z and 
z I—)- e'^'^'^^^^z. A graph of the corresponding Hausdorff dimension 
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6 ^ij^e) = dim//(Aoo(6')) can be found in [T2l Fig. 3]. The fam- 
ily of limit measures {/X6»}o<e<27r/3 represents a continuous transition 
(via singular continuous measures) from the purely atomic measure 
/io := \{5zo + + (^22) to the Lebesgue probability measure /i27r/3 with 
d/i27r/3 := d0/(27r): 

Proposition 2.2. Let 6n,0 E [0,27r/3] and On 0. Then the limit 
measures fig^ converge weakly to fig. 

Proof. The weak convergence in the case 6 = follows from the sym- 
metry of the configuration and the fact that the primitive cells -Dq^^ 
D^f\ D^^ shrink to the points 20,^1,-2^2, as 6n 0. For 6^0, the 
proof of [T2I Thm. 3.1] implies not only the continuity of the Hausdorff 
dimension, but also the weak convergence of the limit measures [131 
Thm. 1.4]. For ^ = 27r/3, Fq is a finitely generated Fuchsian group 
of the first kind. As noted at the end of [15J, the corresponding limit 
measure agrees with the normalized Lebesgue measure. □ 

3. McMullen's algorithm 

McMullen's eigenvalue algorithm [12j provides a very effective nu- 
merical method to compute the Hausdorff dimension of very general 
conformal dynamical systems. For the reader's convenience, we briefly 
recall this algorithm for the special case of a hyperbolic reflection group 
F = Fgg ...^gj.. The algorithm can be extended to obtain arbitrarily good 
weak approximations of the limit measure /ir by atomic measures /ip, 
allowing us to obtain highly accurate values for the moments of fir- 
This high precision is necessary to perform the numerical spectral anal- 
ysis of the CMV matrices in Section O 

3.1. The algorithm. Let e > be given. We start with the prim- 
itive cells D^^^ associated to the geodesies gi and the corresponding 
Euclidean reflections pi in the circles Ci. Next, we generate all cells 
of generation one, {dI^\ d\^\ . . . , D^f^^f^_^^^} , via ps{Df'^) with s 7^ / 
[the 6 second-largest /blue geodesies in Fig. At step n + 1, we start 
with a finite list of cells {dI^\d^\...} and build up a new list of 
cells {D^~^^\ D^^^^\ . . .} by the following procedure. If the radius of 
Dj^^ is less than e, we put it back into the list and relabel it Dp~^^^; 
otherwise we replace it by its k children D^J^~^^\ . . . , D^'^^^^^ (i.e. the k 

cells of generation n + 1 contained in D^^). Since the radius of every 
non-primitive cell is at most rj < 1 times the radius of its parent, we 
eventually obtain a list {D^^\ . . . , D^^^} of cells of generations at most 
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and of radii less than e. We note that every cell in this final list has 
a parent of radius not less than e. 

Let zq, . . . , zm e S*^ be the radial projections of the centres of the 
Euclidean disks d\^\ . . . , -D^/^^ circle. McMuUen introduces 

the notation i \^ j ii Df^"^ C pk{D^^^), where k is the index of the 

primitive cell containing d\^\ and defines the sparse (M + 1) x (M + 1) 
matrix T with the entries 



Let T denote the matrix obtained by raising each entry of T to the 
power (i G M, {T'^)ij = Tf-, and A (A) denote the spectral radius of 
the matrix A. The underlying dynamical system implies that T'^ is a 
primitive non-negative matrix, and we can apply the Perron-Frobenius 
theorem. Note that \{T'^) can be effectively computed via the power 
method. 

We note that the construction (and variables) above depend on the 
choice of e > 0. To emphasize this dependency, we sometimes write, 
e.g., z\'''^ and M(e) instead of zi and M. Let be the unique posi- 
tive number such that A(T'^^) = 1; note that our is denoted an in 
McMullen's article, which here denotes the Verblunsky coefficients in 
Section [i Then ^ Thm. 2.2] 



It is clear from the construction that the sample points Zi and their 
weights wi respect the 6-fold symmetry of Fg, viz., if zi is a sample 
point with weight (i.e. corresponding entry in the Perron-Frobenius 
eigenvector) w;, then so are ^, e^'^'^^^^zi and e^^'^'/^2[. 

3.2. Approximations of the limit measures. Let w'^\ . . . , w^^ , 

with M = M(e), be the entries of the Perron-Frobenius eigenvector 
of normalized so that Y2i '^l'^^ = 1- entries are positive and can 
be considered as approximations of the values /ir(5'^ fl dI^^). In fact, 
we have 

Proposition 3.1. For every e > 0, let 




if i h-^ j, 
otherwise. 



(3) 



lim4 = dimH(Aoo(r)). 



(4) 




(5) 
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be the atomic measure supported on the points zf' (obtained by Mo- 
Mullen's algorithm) with weights w\'^\ Then the probability measures 
/ip^ converge weakly to the limit measure fir as e —> 0. 



Proof. Writing Zi 



and wi 



w 



(^) 



for conciseness, first note that for 



Zg in D 



(0) 



j G {0,1,2}, the corresponding component of the Perron- 



Frobenius eigenvector satisfies 



(6) 



1:1^ 



since = \p'j{,Pj{,z))\ ^. By the uniqueness of the hmit measure 

[Thm. I2.1( iv)]. we only need to show that every weak hmit of yUp^ 
satisfies the transformation property ([T]). It suffices to prove this for 
the generators pi. We discuss the case 7 = po- 

Let 5 > and g G C{S^) be given. Assume first that the intersection 
supp^f n Df^ is non-trivial only for / = 0. Let do = dimj:^(Aoo(r)). 
Since the set Ai = {J^ of'^ has a positive Euclidean distance to both 
the centre and the boundary of the primitive cell D^q \ there exists a 
Cin^ > such that l/Cinv < \p'j{z)\ < Cinv for aU z e Ai. Using (g]), 
there exists > such that, for all z ^ Ai and all e < eg, 

6 



iPo(^)r-iPo(^)i 



do I 



< 



1^1 



(7) 



Since g is uniformly continuous, there exists ei > 0, such that 



\9[z) -g[z 



6 

<2 



V|z - z'\ < 2ei. 



Let e < min{eo,ei} be fixed. We conclude from ([7]) that 



(7(po(^))|p[,(^)|*d/i[f)- / g{z)dfi 



< 



M 

z=o 



9ipoizi))\poizi)\'^'wi ~ / giz)dfi- 



(9) 



Let us assume that all cells -D;^^ contained in the primitive cell D, 



(0) 



are numbered from D^^^ to D)^J ; note that and Mq both depend on 
e. Since g{po{zi)) 7^ only if Z 1— )■ j for some j G {0, . . . , Mq}, because 



only then we have po{zi) G D^^^ C -Do"^ we have 



1(0) 



M 



Y,9{Po{zi))\p'o{zi)\' 



Wl 



Mo 

9Mzi))\p'oizir^wi. (10) 
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Using the fact that /ip^ is an atomic measure and using ([6]), we have 

Mo ^ 

{zi)\'^wi. (11) 

J=U I: l^j 

Putting these together and using (I6l) again for the first inequahty. 

M 



I- Mo 
"^•5' j=Ol:l^j 



f2 9iPoizi))\Poizi)\'''wi - / giz)dfi 



Mo 



Combining ([9]) and ( IT2|) . we obtain 



6 



Mo 



6 



(12) 



< - > Wi < -. 
2 ^ J - 2 



j=0 



I g o po{z)\p',{z)\'^ dp^^ - [ giz)dpP 



< 6. 



(13) 



Now assume that supp (^flZ^o^^ = 0. Then we obtain the same estimate 
( 1T3|) . by applying the above arguments to h{z) = go pq[z)\p'q{z)\'^° and 
using |po(2;)| = \p'q{po{z))\~^ . The general case is obtained by choosing 

a partition of unity Xi; X2 ^ C{S^) with supports disjoint to lj«>i ^ 

and -Dq°\ respectively. 

3.3. Moment computations for SPP. Given a probability measure 
/i on S*^, we have the Hilbert space L^(5'^; p) with inner product 



(0) 



□ 



{f,g):= / f{z)g{z)dp. 



(14) 



Now let 9 G (0, 27r/3) be fixed and consider the limit measure pg of SPP 
constructed above. Associated to pg and L'^{S^;pg) are its moments 



Z ^ dpg. 



(15) 



(9) 



The symmetries of the underlying classical system imply that c 
c^^l G M, and that only every third moment is non-zero [cf. flTK]) below]. 
These moments encapsulate much of the information contained in pg, 
e.g., knowledge of them is sufficient to "quantize" the classical dynamics 
of hyperbolic reflections to the unitary CMV matrices of Section [51 

We obtain approximations for the moments by computing the mo- 
ments of the atomic measure [cf. ([5])] 

Mie) 



(16) 
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for small enough e > 0. With = e^'^^^^ and recalling the 6- fold 
symmetry of the points Zi, we rewrite the points z^^^ as (g-^'^^e^^'^'i- 
with J e {0,1,2}, < / < M/6 = M(e)/6 and 

0<4U...< < ^/2- (17) 
For every fixed /, the Perron-Frobenius entries corresponding to the six 
points (^^^^e^^'^i agree, and we denote their value by m\^^ > 0. Then 
the following holds. 

Proposition 3.2. The moments cl of the limit measure ne of SPP 
satisfy 

= 4?+i = (18) 

and 

M/6 

eg) = (-1)^' lim ^ 6mf ) cos [3j0f . (19) 



Proof. Let Cs denote the sth moment of the measure fig . For brevity, 

i^i"^ and (pi = (l)f 



we write m/ = mf^ and 0/ = The 6- fold symmetry yields 



A//6 

^(e,e) ^ 2m;{cos[(7r/3 + 0;)s] + cos[(7r/3 - ^^s] + (-1)^ cos(s0z)} 




(20) 



if s = ±1 mod 3, 
Y^i 6mi cos(s0/) if s = mod 3. 

By Proposition 13 -H we have fif^ fie, finishing the proof. □ 

Figure [3] shows the numerically computed [6] moments for limit mea- 
sures fie of different opening angles 9. For clarity, we actually plot- 
ted {—l)^cf^, corresponding to the moments of fie rotated by vr [1| 



(0) 

(1.6.63)]. For small 9, the moments c^j appear to oscillate with a 
higher frequency with an amplitude modulated by another lower fre- 
quency. This is not unexpected since, for small 9, the cells of the second 
generation are tiny. Since all points zi fall into these cells, the positions 
of the second generation cells control the general oscillation behaviour 
of the moments. As 9 increases, the cells of every generation grow 
and the distances between successive cells of a given generation de- 
creases, with higher-generation cells contributing different frequencies 
to the graph of the moments. The graph seems to become increasingly 
erratic. 
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Figure 3. Numerically computed non-zero moments 
{-Ifcf^ for 9 = 30°, 50° (top row), 70° and 100° (bot- 
tom row). The horizontal axis is k and e = 10~^°. 

4. Least eigenvalue of hyperbolic surfaces 

4.1. Spectrum of Sda^di,d2 symmetries. Given three lengths 

c/q, c?!, > 0, there is (up to isometries) a unique genus-0 hyperbolic 
surface S^^ of infinite area with three cylindrical ends (funnels) 
whose short geodesies are of lengths 2(io, 2(ii and 2^2 [Fig. |4^]. Refer- 
ring to Fig. Hb, this surface is the oriented double cover of D/Fg^^gj 
where dj is the hyperbolic distance between two geodesies. The three 
geodesies gj, together with the geodesies realising the distances dj, de- 
fine a unique (up to isometries) hyperbolic right-angled hexagon. The 
three heights of this hexagon are concurrent Thm. 2.4.3], and we 
choose their intersection point to be the origin of D. The angle 99 and 
the hyperbolic length t in Fig. HJd can be obtained using trigonometric 
identities for trirectangles and pentagons [H Ch. 2], giving 

COS if 



^51 + 51 + 25^5^52^51 + 51 + 25o5^52 
tanht = 5^52 



51 + 51 + 51 + 25o5i52 



(21) 



[51 + 5151 + 25o5i52){5l + 51 + 25^5i52) 
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where 5j = coshrfj. The Euchdean distance from the origin of 3 is 
X = tanh(t/2). These resuhs allow us to compute explicitly, for given 
distances do,di,d2, the three primitive cells needed for McMullen's al- 
gorithm. 



Figure 4. (a) The hyperbolic surface Sdo,dj^,d2- (b) The 
corresponding hexagon with concurrent heights. 

Let D{do,di,d2) := Aoo(X g^^g-^^g^) denote the Hausdorff dimension of 
the limit set, and let \o{SdQ,di,d2) be the bottom of the spectrum of 
the (positive) Laplace operator on the surface Sdo,di,d2- The following 
result is well known and relates the above two functions. (The result 
holds generally for geometrically finite hyperbolic manifolds of infinite 
volume.) 

Theorem 4.1. f2R Thm. 2.21] We have, with the notions from above, 



and Xo{Si^^ij^^i^) < 1/4 is the eigenvalue of a positive L"^ -eig en function 
if and only if D{dQ,di,d2) > 1/2. 

The spectral interpretation of the Hausdorff dimension was extended 
by Patterson in [16] (see also [21 Thm. 14.15]) to the case D = D{do, di, ^2) < 
1/2, in which case D represents the location of the first resonance. 

Fixing a total distance c? > 0, we consider the function {do, di, ^2) ^ 
D{dQ,di,d2) with do + di + d2 = d. One might think that the point 
Po = {d/3, d/3, d/3) would be a unique global minimum of the Hausdorff 
dimension subject to this total distance restriction (in fact, this was 




(a) 



(b) 




(22) 
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our initial guess). Numerical computations showed, however, that the 
function D{do,di,d2) is more comphcated and has, in fact, precisely 
four global minima at Pq and at 3 other points, 

/2d d d\ /d 2d d\ a ^ '^^\ 

iT'e'eJ' U'T'eJ U'e'Tl 

This observation lead to the following rigorous result: 

Proposition 4.2. Let a,b> 0. Then we have 

D{2b,a,a) = D{2a,b,b). (23) 

We note that Proposition 14.21 shows that, for (i > 0, the function 
t I— >■ D(t, ^) is symmetric with respect to d/2. In Fig. [H this 
means that D is symmetric about the midpoint along the line AB. 

Proof. The hexagon with the three geodesies go,gi,g2 with distances 
a, a, 2b is shown in Fig. [5b. The configuration is obviously symmetric 
with respect to the horizontal geodesic g^. Let pj denote the hyper- 
bolic reflection in gj and Dj be the closed disk with gj = dDj. Then 
^90 91 92 = (Po,Pi,P2) =: To, and a fundamental domain of Fq is given 
by' ' 

J^o = {zeC: Re{z) > 0}\(Do U D^). (24) 
The geodesic g^ in Fig. 03 is the reflection of g2 in gi. Therefore, 
we have for the corresponding reflection = piP2Pi- A fundamental 
domain of F = (po, Pi, P2, Ps, pi) is 

F ={zeC: Re{z) > 0, lm{z) < 0}\D2. (25) 

Since J-q = J-" U pz{J^) we have [F : Fq] = 2. The geodesies g2-,g^,gi 
in Fig. [5]d have distances 6, 6, 2a. Completely analogously, noting that 

^92,53,94 = {P2,P3,P4) =■ Ti has 

Ti = {zeC: lm{z) < 0}\(L>2 U D4) (26) 

as a fundamental domain, we have po = P3P2P3, and J^i = U pi{J^). 
This implies, again, [F : Fi] = 2 and, therefore, 

dim^^Aoo(Fo) = dimHAoo(ri) = dimHAoo(f ). (27) 

□ 

Remark. The spectrum of a geometrically finite hyperbolic surface of 
infinite area consists of an absolutely continuous part [1/4, 00) without 
embedded eigenvalues and finitely many eigenvalues of finite multiplic- 
ity in the interval (0, 1/4); see jSHII] or ^ Ch. 7]. 
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Assuming D{2a, b, b) > 1/2, we conclude from Thm. I^TTl and Prop. 112] 
that Xo{S2a,b,b) = Ao(5'26,a,a), i-G. S2a,b,b and S2b,a,a havc the same lowest 
eigenvalue. It is natural to ask whether these two surfaces are iso- 
spectral. We obtain S2b,a,a from S2a,b,b by first cutting S2a,b,b along the 
geodesies a, (3 in Fig. [S^i, unfolding it onto D\(i5o U D2 U -D4 U D^) 
=: Ds in Fig. [5]d, and then gluing the boundary geodesies g^^g^ and 
g2, (74 together. Let pj be the reflection along geodesic gj as before, with 
pi,P3 also considered to act isometrically on S2a,b,b and S2b,a,a- Since 
pi, p3 and the Laplacian A commute, we consider simultaneous eigen- 
functions of these operators on S2a,b,b- Let / be an L^-eigenfunction in 
S2a,b,b with eigenvalue A < 1/4 which is even under pi and pa, 

Af = Xf with /opi = / = /op3. (28) 

(At least one such / exists — with A = Aq.) Now any such / G C°°{S2a,b,b) 
can be transplanted to / G C°°{Ds)- Thanks to the symmetry fl28b). 
/ along with all its derivatives agree along the corresponding points on 
go, g5 and g2, g^, so / can in turn be transplanted to an L^-eigenfunction 
/ G C°°(5'26,a,a) with the same eigenvalue. A/ = A/ in S2b,a,a- One can 
carry out the same argument when / o p^ = — / = / o pg in ( 128|) . as- 
suming such an / exists. This argument shows that some eigenvalues 
of these two surfaces coincide. 
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4.2. Example Sdo,d-^,d2 with do + di + d2 = 3. Identifying the triple 
{do, di, d2) with the point {do + die^'^'/^ + rfse^'^'/^) i/y/?> e ]R^ we rep- 
resent the domain of D{do,di,d2) by an equilateral triangle T C 
centred at with heights see Fig.[TJ In Fig.Elwe plot D{do, di, ^2) 

restricted to do + di + d2 = 3, computed using McMuUen's algorithm. 
When any dj is small (i.e. for points close to dT), it is difficult to 
compute D numerically using McMuUen's algorithm even with 256-bit 
arithmetic; we have therefore left out the blank regions near dT (thick 
black triangle on the plot). 




We found the minimum of D to be D{1, 1, 1) = 0.56996 ... > 1/2, 
which as noted above is also attained at three other points (2, 1/2, 1/2), 
. . .. We remark that D is very flat for much of T (red region in the plot), 
only increasing rapidly near dT (where our numerical computation 
breaks down). The four global minima of D can also be seen in Fig. [TJ 
which shows the bottom of D at an expanded vertical scale. In view of 
the apparent smoothness of the function D in Fig. [TJ it is natural to ask 
whether there are explicit formulas for the gradient or the Hessian of D. 
From the shape of the graph (the part reliably computed), we believe 
that D attains its global maximum at the vertices and the midpoints 
of dT. 
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Figure 7. The bottom part of D{do,di,d2); note scale. 



The value of D at these points can be computed using (123|) and 
the results of Baragar [1], who described an algorithm for the Haus- 
dorff dimension of reflection groups Tg^^g-^ g^ with two distances equal 
to zero. Our {do, di, ^2) = (3, 0, 0) corresponds to a = {l + x)/(l — x) = 
2.35240... in P with x = tanh2(3/4). Using [H Table 1], we ob- 
tain 0.225 < Ao(S3,o,o) < 0.226 and consequently 0.654 < £'(3,0,0) = 
D(0,3/2,3/2) < 0.659, which we believe bounds the global maximum 
of D in r. 

Remark. Our numerical computations indicate that the situation does 
not change qualitatively when we choose a much larger total distance 
d (we also carried out detailed calculations for the case d = 6 with 
the global minimum 0.334541 . . ., and obtained graphs very similar to 
Figs. [6] and [7]). For any e G (0, eo), there exists a total distance d = d{e) 
such the four global minima of D under the restriction do + di + d2 = d 
agree with e. Again, for a large part of T, the graph of D should be 
relatively fiat. 

It seems plausible and there is strong numerical evidence (although 
we do not have a proof) that one has 

D{do, di, d2) < D{d'Q, d[, 4) if ^ < dj for j = 0, 1, 2. (29) 

Assuming continuity of D up to the boundary dT, the graph must be- 
come very steep close to the boundary, because McMullen's asymptotic 
limc^.oo D{c, c, 0) = 1/2 [12j Thm. 3.6], together with the monotonicity 
property fl2^ . imply that D assumes values > 1/2 at all boundary 
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points dT: For (a, h, 0) G dT with a > 6 > we have 

D{a,h,Q) > D{a,a,0) > hm D{c,c,0) = 1/2. 

c— >oo 

5. OPUC AND CMV MATRICES 



(30) 



5.1. Overview. We review briefly the properties of OPUCs and CMV 
matrices as relevant to the current work, referring the reader to [T8142T] 
for more details. As before, D C C denotes the open unit disk and 
= (9D. Given a probability measure /i on 5*^, supported on an infinite 
set, and the inner product (HM . let /^)}fcLo be a set of monic 

polynomials orthogonal with respect to fi (with the usual convention 
that $fc(-2; yw) = z'' + lower-order terms); for brevity, we will often write 
$fc(z) for ^k{z] fi) when there is no confusion. As do real orthogonal 
polynomials, {$fc} satisfy a recurrence relation 

^k+iiz) = z^kiz) - ak^iz) (31) 

called Szego's recursion, where the Verblunsky coefficients ak can be 
shown to lie in 3. The reversed polynomial is 

n{z) := z'^^ilM (32) 
or, with $fc(z) = J2'j=o bjZ^, 



Ej=0 (^k-jZ^ 



This implies that $^(0) = 1 for all k and, together with ( !3T|) . 

ak = -$fc+i(0). 



(33) 



(34) 

fc T oo 
A:=0- 



In practice, one can compute {^k}kLo using Gram-Schmidt on {z } 
We note that {$fc} may or may not form a basis for L?'{S^\ /i); see [20 
Thm. 2.2]. 

If, on the other hand, we apply Gram-Schmidt to {1 
. . .}, we get orthonormal polynomials {xo{z) , xi{z) , X2{,z) , . . .} which 
do form a basis for L'^{S^;fi). The CMV matrix associated to the 
measure n is the matrix representation of the operator f{z) — )■ zf{z) 
on L'^{S^; /x). It has the semi- infinite pentadiagonal form 



C 



( ao 
Po 




V : 



aiPo 
— aiOo 

P2P1 




PiPo 
-Pitto 








a3P2 
-a3«2 
a4P3 






P3P2 
-P3a2 



(35) 
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where pk = a/1 — \otk\'^- We note that Jacobi matrices, obtained in a 
similar way for orthogonal polynomials on the real line, are tridiagonal 
matrices. As in the case of orthogonal polynomials on the real line, an 
important connection between CMV matrices and monic orthogonal 
polynomials is 

<l>„(2) = det(z/-C(")) (36) 

where C*-"^ is the upper left n x n corner of C. 

If |q;„_i| = 1, C decouples between {n — 1) and n as C = C*^"^ © 
C, where the upper left corner is an x n unitary matrix C^"-* = 
C*^"'^(ao, «i, • • • , «n-i) and the remaining block C{an, On+i, • • .) is a (semi- 
infinite) CMV matrix. This suggests that a unitary n x n truncation 
of a CMV matrix can be obtained by replacing a„ G D by /3 G d3. 
The truncated CMV matrix has as characteristic polynomial 

<^n{z;(3) = z<!>r.-i{z)-f3K-iiz) (37) 

whose zeros are all simple and lie on dH) = S*^, a fact which will be 
convenient below. The polynomials $„(2;;/3) are called paraorthogonal 
polynomials. 

Given a set of moments Cj^lfi), we can recover its generating measure 
H in the classical limit as follows [HI Thm. 2.2.12]. Given /3 and k, let 
{zj}'j^i C (9© be the zeros of ^k{z; (3) and define the atomic measure 

k 

Cf := y 5^^ (38) 

where ^Pi{z) = $j(2;)/||$j||, and || ■ || denotes the norm in L'^{S^\p). 
For any choice of /32, /Ss, • • • G 9D, C,^^^ converges weakly to /i. 

The weight |y9j(zj)p) ^ is known as the Christoff el function [HI 
p. 117ff]. 

5.2. OPUCs for SPP. Let us now return to McMullen's SPP. Using 
the moments c[.^^ computed in Section 13. 3[ we used Gram-Schmidt to 
construct the orthogonal polynomials $fc(z) = (^k{,z]ii0) associated to 
the measure /x^, which for convenience we rotate by vr. Henceforth, 
by po we mean this rotated measure. Since all moments are real and 
'^sfc-i — ^sfc+i — 0' polynomials $3^(2;) are, in fact, polynomials in z 
with real coefficients, and $3^+1 (z) = z(^^k{,z) and (^^k+2{,z) = z'^^^kiz). 
It then follows from that the Verblunsky coefficients are all real 
with a^k = «3fe+i = and 03^+2 G (—1, 1); this also follows from the 
symmetries of /i and [HI (1.6.66)]. 
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It is convenient to introduce monic polynomials qk{z^) '■= ^3k{z) 
with Verblunsky coefficients 7^ := a3k+2, in terms of which (1^ reads 

qk+i{z) = zqk{z) - -fkq*k{z). (39) 

In analogy with (1371) . we define the paraorthogonal 

qk{z; /3) = zqk-i{z) - (3ql_,{z), (40) 

where we choose now /3 = ±1 depending on the sign of 7^-1. Since 
qk{z^] P) = ^?,k{.z] (3), the zeros of qk{z; (3) all lie on 9D. Moreover, since 
qk has real coefficients, its (non-real) zeros occur in complex conjugate 
pairs. For our numerical computations, we used qk exclusively in place 
of $3fc. Since the OPUCs depend on the underlying measure /z only 
through the moments Cj, it is clear from the definition (fT5|l of the latter 
that {qk}'k=o the OPUCs for a measure jle where 

d^e(0) = 3d^e(0/3) (41) 

for <p e [-7r/3,7r/3). 

Having computed qk and 7^ = — q'fc+i(0), the accuracy of the numer- 
ical computation can be checked using flHT]) . by ensuring that the error 
(which is zero for exact computation) 

'■= J2j \Qk+l,j — qk,j-l + 'ykqk,k-j\, 

with qkj the jth coefficient of qk (with qk,-i = 0), remains small. We 
found that high-precision arithmetic, both for the moments Ck and 
the subsequent computations involving q^, are crucial to control the 
error. For 9 = 119° and e = 10^'', computations using a 256-bit "quad- 
double" precision gives us /c ~ 280. 

5.3. Observations. We now present a few numerical observations on 
the spectral properties of CMV matrices. In this section, we work 
exclusively with the symmetry-reduced polynomials qk and the measure 
fig introduced above. 

First, for every 6 G (0, 27r/3), the Verblunsky coefficients are all 
negative, 7^ < 0. Seen in the light of the formula for Verblunsky 
coefficients for rotated measures [IHl (1.6.66)], for any measure /i, 

7.(/i("))=e-('=+^)7fe(^) (42) 

where d/i*^"-* (</>) = dfi^cp — a), our observation means that yUg belongs to 
a family of measures whose Verblunsky coefficients are, possibly after 
rotation, all negative. To obtain paraorthogonal polynomials qk{z; 13) 
corresponding to unitary truncations of the CMV matrix, it is therefore 
natural to take (3 = —1 (but see the effect of the choice of /3 at the end 
of this section). 
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Figure 8. Monotonicity of Verblunsky coefficients: 
Plot of log(l + 7fc) against 6 for (from top) k = 1,2, 4. 

Secondly, for every k fixed, 'jkif^e) increases monotonically from — 1 
to as ^ goes from to 27r/3. Figure [H] plots log(l + 7fc) for k = 1, 2 
and 4 against 9, showing the rapid convergence of 7^ to — 1 as 6^ — > 0. 

In the course of our numerical computations, we found that small 
opening angles 6 only allow us to obtain a few polynomials qk reliably, 
while larger 6 allows us to obtain more polynomials. Since the numer- 
ical instability of the Gram-Schmidt orthogonalisation of the polyno- 
mials l,z,z'^, ... is closely related to the numerical singularity of the 
Toeplitz matrices 

= {Cjli)o<i,j<n = ((^'', ^*)e)o<ij<n, (43) 

we expect that the determinant Dn{dng) := det T^"'"''^'' decreases mono- 
tonically to as 6' — )■ 0. In fact, this would follow from the monotonicity 
of the individual Verblunsky coefficients by the identity [HI §1.3.2 and 
(2.1.1)] 

n~l 

D^{df,g) = l[il-\^,{fig)\'r~^. (44) 
j=0 

In Figure [H we present our numerical computations for 6 = 119° 
with e = 10"'^ (size of the largest cell). The thin (green) line is the 
distribution function of the classical measure /iii9([0, 0]). As noted 
earlier, dfig tends to the Lebesgue measure d0/(27r) as 6* — )■ 120°, but 
significant "gaps" (i.e. intervals I with = 0) remain even for 6 = 
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Figure 9. Solid (green) line: (classical) measure /ing, 
red dots (overlapping the line): C,^'^^^\ blue dots: idz*-^^^''. 

119°. Let Zj be the zeros of the paraorthogonal 5'28o(-2) = ordered 
from = to 27r. The blue dots in Fig. |9]plot the "integrated density 
of zeros" measure 

while the red dots plot the "Christoffel-corrected" measure 

^(280) 

fined in fl551) . One can see that (^(2*^°) overlaps the classical measure /ing 
to within plotting line thickness. This agreement confirms the accuracy 
of our numerical computations. We note that one may find zeros in 
the spectral gaps S'^\supp/ie, but with smaller probability as — > oo 
and the Christoffel function at these "spurious" zeros is small. 

The idz plot also suggests that in the limit /c — )■ oo, zeros appear to 
accumulate at the gap edges, most visibly for the main gap centred at 
= 27r and the secondary gap starting at vr. As noted above, there 
could be zeros inside the gaps; here they are most visible near 27r/3 
and 47r/3 [see also Fig. [10] below]. However, their contribution to the 
measure C,'^^'^ is weighted down by the Christoffel function. 

In Figure [TOl we plot, for 9 = 119° and /3 e S^, the migration of 
the zeros of g4o(-2; /?) as arg/3 varies from to 2n, superimposed on 
the support of the classical measure fi^g (since supp fi^g has Lebesgue 
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TT arg Zk 

Figure 10. Migration of zeros of 540(2;; /3) as a function 
of (3 for 9 = 119°. Horizontal axis is aigZk and the 
vertical axis is arg /3. Each vertical "line" is in fact made 
up of 500 disjoint dots. Faint (green) background is the 
approximate support of the classical measure /^ng. 

measure zero, the "solid" background arises only from plotting limi- 
tations). This plot confirms Thm. 1.1 in [22], which states that any 
interval [a, b] ^ supp /x, represented by a white background, contains at 
most one zero of q{z;fi,(3). It also confirms Thm. 1.3 in [22], stating 
that for any /3, f3' G S^, the zeros of qk{z] (3) and of qk{z; (3') strictly 

interlace. The graph also shows that each zero z^f'' is monotone in 

arg /3 and travels from z^^^ to z^^^^ as arg /3 increases from to 27r. 
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